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Abstract 

In critical as well as in non-critical string theory the partition function reduces 
to an integral over moduli space after integration over matter fields. For non- 
critical string theory this moduli integrand is known for genus one surfaces. The 
formalism of dynamical triangulations provides us with a regularization of non- 
critical string theory. We show how to assign in a simple and geometrical way a 
moduli parameter to each triangulation. After integrating over possible matter 
fields we can thus construct the moduli integrand. We show numerically for c = 
and c = —2 non-critical strings that the moduli integrand converges to the known 
continuum expression when the number of triangles goes to infinity. 
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1 Introduction 



Non-critical string theory, or two-dimensional Euclidean quantum gravity coupled 
to matter, has been a fruitful laboratory for studying aspects of string theory as 
well as quantum gravity. Solving these theories one has had the advantage to have 
both a lattice version of the theory and a continuum field theory formulation. 
The lattice version has been denoted the dynamical triangulation model (DT) 
or the matrix model of non-critical string theory. It can be solved, basically by 
combinatorial methods, exemplified by the use of matrix models. The explicit 
solution of the continuum model uses the fact that the Liouville theory is a 
(special) conformal theory. For observables which can be calculated by both 
approaches agreement is found. 

Since non-critical string theory is described by conformal field theory we know 
that conformal invariance is implemented. However, the precise manifestation 
of this invariance, and how it is related to the moduli space of the underlying 
surfaces (i.e. to the part of the surface geometry left invariant under conformal 
transformations) has been rather limited. That dynamical triangulations or ma- 
trix models contain precise information about the moduli space is on the other 
hand obvious. For instance it was shown in [1] that the resolvent of the matrix 
model has an expansion in terms of so-called moments, where the coefficients in 
the double scaling for the genus h terms are precisely the intersection indices of 
moduli space for genus h Riemann surfaces with any number of punctures. It was 
also shown how these matrix models in the double scaling limit can be related to 
the Kontsevich matrix model which provides a representation of the generating 
function for these intersection indices [1, 2]. Much work has later expanded on 
these results [3], but we are in these approaches far from the naive and simple 
idea which was the starting point for DT and matrix models, namely that DT 
provides a regularization of the path integral of non-critical strings and thus the 
moduli parameters should appear in the same simple way as they formally appear 
in the string path integral. 

One problem when using matrix models is that most observables which can 
be calculated analytically are of global nature: integrated correlation functions of 
matter fields or the partition function as a function of various boundary states. 
Further, the comparison with continuum results is only possible for the simplest 
topologies of the surfaces: spherical topology, disc topology and (only recently [4]) 
cylindrical topology. It is simply difficult to perform calculations in the frame- 
work of Liouville theory for higher genus surfaces. However, there is a narrow 
window where one can test in more detail if the matrix models, or the framework 
of dynamical triangulations, actually agree with the continuum expressions for 
higher genus surfaces in the naive way mentioned above. We know the partition 
function for genus one surfaces of non-critical strings expressed as an integral over 
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the moduli parameter of the torus [5]. Integrating out the matter fields, using 
the conformal anomaly, worldsheet conformal invariance ensures that the remain- 
ing integrand of the partition function is only an (explicitly known) function of 
the moduli parameters. We can now in principle compare this integrand to the 
integrand which we obtain using dynamical triangulations as a regularization of 
the non-critical string theory. Which integrand do we obtain using dynamical 
triangulations? For each triangulation we have a continuous, piecewise linear ge- 
ometry. To each such geometry we can associate moduli parameters (as described 
below). Each geometry will appear with a certain weight which depends on the 
matter we have coupled to the surface. If we have N triangles we will in this 
way obtain a number of points in moduli space and when goes to infinity we 
expect these discrete points converge to a density distribution which should be 
proportional to the integrand of non-critical string theory. The possibility of per- 
forming such a comparison was pioneered by Kawai and collaborators [6]. They 
found good qualitative agreement between the non-critical string integrand and 
the density constructed from dynamical triangulations. The purpose of this arti- 
cle is to improve this test, making it quantitative, and also present a more general 
setup than the one used in [6]. We also show how one can measure the moduli 
parameters for higher genus triangulations, but unfortunately there is presently 
no theoretical calculation with which we could compare such results which would 
be easy to generate numerically. 

2 The setup 

In string theory we can integrate out the conformal matter fields on the worldsheet 
by the conformal anomaly. This is true also for non-critical string theory. The 
partition function of e.g. a bosonic string in d dimensions can be written as 

Z^'^) = J dllin) J V^(f)V^bV^cV^X^ ^-S{X,g)-S{b,c,9)-Sj.{^,9) ^ 

where the integration over worldsheet geometries of genus h is implemented by 
integrating over metrics gafi, gauge fixed to conformal gauge ga/s = e-'^gajsiTi)-! 
ga/3{'Ti) denoting a fiducial background metric depending on (for h > 1) 3/i — 3 
complex moduli parameters Tj parametrizing the complex structure of the genus 
h Ricmann surface and d/x(rj) being a modular invariant measure. The b, c denote 
the Faddccv-Popov ghosts associated with the partial gauge fixing to conformal 
gauge, (p denotes the Liouville field and finally SL{4>,g) is the Liouville action: 

SL{^rg) = ^j d''i^)({do,^f + QR^ + ii^^^) (2) 
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The constants in the Liouville action are chosen such that the conformal in- 
variance is maintained in the quantum theory in accordance with the bootstrap 
approach of David, Distler and Kawai [7], i.e. 

Q = V(25 - d)/6, g = + p. (3) 

For a fixed value of the moduh parameters one can in principle perform the 
integration over the bosonic fields and the ghost fields. One obtains contributions 
from zero modes and appropriate determinants. However, only in the case of 
genus one is the integration over the zero mode trivial, the reason being that 
/ d^^^/^R = 0. In this case one obtains just the standard result from critical 
string theory [5]. Thus, if we consider the situation where we fix the world sheet 
area to be A, i.e. trade the cosmological constant jj, in (2) for the area A, we 
obtain 

Z^^='\A) A-' [ ^F{Tf-\ (4) 

In (4) c denotes the central charge of the conformal matter theory coupled to 2d 
gravity, r = ri + iT2 the complex moduli parameter of the complex structure we 
are integrating over and M. is the region in the complex plane to which r belongs. 
Finally F{t) is given by: 

oo 

F(t) = Ts-'/'e'^^/^ n 1^ ~ e''^^"^|-'. (5) 

n=l 

We want to test how well the non-critical string theory regularized by DT 
approximates this formula. We will further restrict ourselves to the cases c = 
and c = —2. We do so in order to obtain the best statistics in the numerical tests 
we perform. The case c = is chosen to avoid to have to update matter fields in 
addition to the geometry, when we use Markov chain Monte Carlo simulations to 
generate configurations. The case c = — 2 is chosen in order to have a non-trivial 
matter system coupled to gravity which can be easily handled numerically. For 
c = —2 there exists a recursive algorithm which allows us to avoid construct- 
ing Markov chains and generate directly random triangulations with the correct 
weight, including the matter fields [8]. 

The set-up is thus the following: Using DT we can fix the area, hke in (4). 
This is done by fixing the number of triangles, since all triangles are viewed as 
identical in the DT formalism. A triangulation can thus in principle be viewed 
as representing a piecewise linear geometry. As we will show below there is a 
natural way to associate with such a piecewise linear geometry a moduli parame- 
ter. For the ensemble of piecewise linear geometries constructed from N identical 
equilateral triangles, each triangulation with the topology of a torus, we will then 
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have an approximation to the continuum distribution r2~^-F(r)'^~^. We will in- 
vestigate if and how the corresponding distribution of the moduli parameter r 
constructed from the DT ensemble characterized by N triangles converges to the 
continuum distribution ior N ^ oo. A priori it is not clear that there should 
be a convergence at all. Given a continuum area A, we usually think of the lat- 
tice area as A = Na^^/^i/A, a being the link length of the triangulation. Thus 
keeping A fixed and taking N ^ oo dictates how one should think of the lattice 
link length a — ?> 0. However, the corresponding piecewise linear surface will in 
general not be anything like a smooth surface. In fact we know now that it with 
probability one will be fractal with a Hausdorff dimension different from two [9] . 
The situation is thus quite similar to the one encountered for the ordinary path 
integral of a particle: with probability one a path in the path integral is nowhere 
differentiable and has a Hausdorff dimension different from one. As long as we 
think of the DT ensemble of piecewise linear surfaces in this way one would not 
be too worried about the highly fractal structure of the geometry associated with 
a generic triangulation in the ensemble for N ^ oo, and we know that for global 
observables one obtains identical results using continuum Liouville theory and 
the DT ensemble, taking iV — )■ cxd. However, it is less clear if an assignment of a 
moduli parameter as we are going to do it, by analogy between differential forms 
and "discrete differential forms" on simplicial complexes, even if very natural, will 
work for the "wild" generic triangulations we meet in the DT ensemble. Recall 
for instance that the spectral properties of the Laplacian defined even on the sim- 
plest piecewise linear surfaces are different from the generic spectral properties of 
the Laplacian defined on a smooth geometry, because one should, when analyz- 
ing the spectrum, view the piecewise linear geometry as a geometry with conical 
singularities at the vertices [11]. Since we are actually using properties of the 
Laplacian to determine the moduli parameter, one could be worried whether the 
assignment using discrete differential forms makes sense in the present context. 
However, we will find that the assignment works beautifully! 

Let us start by showing how to determine the moduli parameter r if given a 
two-dimensional Riemannian geometry of a genus one surface. Let d denote the 
exterior derivative and 6 its adjoint (the co-differential) with respect to the stan- 
dard inner product defined on p-forms, p = 0, 1, 2 for two-dimensional surfaces: 



In components the Hodge dual of a p-form i/jp on a n-dimensional Riemannian manifold is 




(6) 



where *ipp is the Hodge duaP 2 
1-forms: 



—p form of ijjp. In particular we have for two 




(7) 
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where 0i = ipaix)dx°' and ipi — %l)a{x)dx°', a = 1, 2 in some coordinate system x^. 
The Hodge Laplacian 

A = d5 + 5d (8) 

maps p-forms to p-forms, and for 1-forms on a genus h surface the kernel is 2h- 
dimensional, constituting the so-caUed harmonic differentials. For a genus one 
surface the vector space of harmonic 1-forms is thus two-dimensional. We can 
find it by solving 

d(t)i = 0, 5(t>i = 0, (9) 

which express that a harmonic 1-form has a vanishing curl and divergence. 

Let us choose two closed curves 71 and 72 on the torus that generate the 
fundamental group. It is then possible to choose "dual" harmonic 1-forms and 
a^, i.e. harmonic forms which satisfy: 

/ o^-^l (10) 

For the torus the moduli parameter can be defined in the following way: by 
the uniformization theorem any metric on a surface is conformally related to a 
constant-curvature metric. For the torus it is a zero-curvature metric, and starting 
from the Euclidean plane we obtain such a flat torus by identifying opposite sides 
in a parallelogram. Let the two sides of the parallelogram be vectors uji and 
U2, represented as complex numbers. The moduli parameter r is defined as 
T = UJ2/UJ1. It is clearly invariant under translations, rotations and global scaling 
of the parallelogram and we can assume that T2 > 0. However, there are finite 
diffeomorphisms of the torus which cannot be obtained continuously from the 
identity, and which we will insist do not change the partition function. These 
transformations do not leave r invariant but will change it according to 

ar + b a.b.c.deZ 

^ h \] u 1 ■ (11) 

CT + d ad-bc^l ^ ^ 

The transformations (11) constitute a transformation group called the modular 
group r and it is isomorphic to SL{2, Z)/Z2. In the integral (4) we are instructed 
only to integrate over points r which cannot be identified by the action of the 
modular group via (11), i.e. Al is a so-called fundamental domain of F in the 



defined as: 



pi 

^'^P = M £ai...a^,6i...6„_^V'"^-"'' dx^' A • • • A dx^^^ 

p\[n — py. ^ ^ 
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upper half-plane: no pair of points within it can be connected by a modular 
transformation and any point outside it can be reached from a unique point inside 
by some modular transformation. A standard choice of fundamental domain Ai 
is: 

^ ■( f ^2 > 0. - i < Ti < and |r| > 1 

T e if < „ 2 ^ . (12) 

I r2 > 0, < Ti < i and |r| > 1 ^ ^ 

We will use this fundamental domain in the rest of the article. Of course it is 
not unique: for example any modular transformation will transform it to another 
fundamental domain. 

Given a metric gab and the curves 7^, r is determined by the formula 

To prove this formula we note that we can choose periodic coordinates and 
x^, both with period 1, such that 71 and 72 run around in the x^ and x^ direction, 
respectively and such that metric takes the form: 

ds^ = e''^''^gabdx''dx\ (14) 

where e'^*^^^ is a local scale factor and g represents a flat metric with determinant 
1. The parallelogram can now be represented in the Euclidean plane as shown in 
Fig. 4 and the g is thus explicitly given by 

the factor I/T2 enforcing det gah— 1- Since harmonicity of the 1-forms is preserved 
under conformal transformations, wc clearly have that = dx^ and — dx^ 
and the inner product {a^\a^) can be read off from the definition (7): 

{a'\a') = g», i.e. {a^\a^){a^\a^) - {a^\aY = 1 

Eq. (13) follows from (15) and (16). 

The T determined by formula (13) lies in the upper half-plane, but not nec- 
essarily in the fundamental domain M. defined in (12). We use (11) to find the 
modular transformation which maps it to A^. 

The formula (13) can be applied directly to the piecewise linear geometries 
encountered in the DT formalism since the notion of p-forms and an exterior 
derivative acting on these forms can be defined in a natural way. For a detailed 
definition for piecewise linear geometries, defined by triangulations where the 
lengths iij of the links between vertices i and j are given, we refer to standard 

7 



(13) 



(15) 



(16) 



Figure 1: The embedding of a triangulation of the torus into the Euchdean plane. 
This particular triangulation with = 280 triangles is taken from the ensemble 
DT{3) and has moduli parameter r = 0.09 + 0.7i. 



textbooks [12]. In the case of DT the formulas become very simple because 
all link lengths are identical, and similarly all areas of triangles are identical. 
Thus for the DT piecewise linear surfaces the definitions become identical to the 
"abstract" definition of discrete differentials used on simplicial complexes to study 
the discrete versions of cohomology as a more general setup than the De Rham 
cohomology of ordinary differential forms. A discrete p-form 0p is in the simplicial 
complex context defined by assigning a real number 0p(cTp) to each oriented p- 
simplex dp. Changing orientation of the simplex changes the sign of the assigned 
real number. The linear space Qp of p-forms has dimension equal to the number 
of p-simplices in the triangulation. For the scalar product corresponding to (6) 
we simply choose 



where the sum is over all p-subsimplices in cip+i and the sign factor signifies the 
orientation of the subsimplex relative to the orientation of o"p+i. Again S, the co- 
differential, is defined as the adjoint of d with respect to the scalar product (17). 
To be explicit, consider the 1-forms. The 1-forms take values on the oriented 
links. Denote the vertices in the triangulation i and the links (ij) when vertices 
i and vertices j are connected by a link oriented from i to j. Similarly triangles 




(17) 



The exterior derivative is defined as 
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are denoted (ijk). Thus we have = ~<P{ji) a-nd 

{d(p){ijk) = + 0(ifc) + 0(fci), = ^ (19) 

where the summation is over all links incident on vertex i. If one wants to 
relate 1-forms defined this way to continuum notation, one can heuristically think 
of the 1-form as 

0(i,)= r dx'^Mx). (20) 



Using (20) and appeahng to Stokes theorem {d4>)(ijk) will be proportional to the 
curl associated with triangle (ijk), and similarly Gauss' theorem makes {S(f))i 
proportional to the divergence assigned to the cell dual to vertex i. We note that 
using the definition (19) it is easy to show [12] that for a triangulation of genus 
h the equations {d(j))(^ijk) — and {S(f))i ~ have 2h solutions, in agreement with 
the continuum result referred to above. 

Finally we define a path 7 in the triangulation as a path of links. Thus we 
have a natural integration of 1-forms along a path: 



E '^fe)- (21) 

{ij)e-y 



We now have all the ingredients for the construction of the moduli parameter r 
for a piecewise linear triangulation belonging to the DT ensemble of genus 1. First 
identify two non-contractible loops 71 and 72 which generate the fundamental 
group. Next solve (d0)(yfe) = and {S(j))i = to find the two-dimensional space 
of harmonic 1-forms. Finally select from this space the two 1-forms and 
which satisfy the discrete analogy of (10), i.e. 

E = (22) 

and use (13) (and (17)) to calculate r. 

Let us finally note that we can use the harmonic 1-forms and to obtain 
an explicit embedding of a DT triangulation into a parallelogram in the Eu- 
clidean plane. In the continuum, choosing coordinates such that the generating 
curves 71 and 72 run around in the and x"^ directions, respectively, we have 
a* = dx*. Thus, knowing a* we can recover the coordinates x^ by integration. 
We now apply this to our triangulation using the two harmonic 1-forms a^j^.^^ we 
have constructed for the triangulation. For each link [kl) we view the vector 
{oiQ^iy Oi^^j^^s^ as representing the link {kl) in the Euclidean plane precisely in the 
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Figure 2: Embeddings of two large triangulations of the torus into the Euchdean 
plane. The triangulations have 150000 triangles. The triangulation to the left 
is picked from the ensemble of c = triangulations (2d gravity without matter 
fields), while the triangulation to the right is from the c = —2 ensemble (to be 
described in Sec. 4). Shown in red are the shortest non-contractible loops. 

heuristic sense of eq. (20). Starting from an arbitrary vertex we can now recon- 
struct the triangulation in the Euclidean plane. Fig. 1 shows an example of such 
a reconstruction. Recall that a function is harmonic if and only if it has the mean 
value property: the value of a harmonic function at the center of a circle is equal 
to the integral of the function along the circle divided by the length of the circle. 
The discrete map from the triangulation to the Euclidean plane is harmonic in 
the sense that any vertex is located at the center of mass of its neighbors. 

Fig. 2 shows other embeddings of triangulations of the torus in the Euclidean 
plane. These triangulations are much larger than the one shown in Fig. 1, each 
consisting of = 150000 triangles. The triangulation to the left is from the c = 
ensemble of triangulations considered in the next section (2d gravity without 
matter fields), while the triangulation to the right is picked from the ensemble of 
triangulations describing conformal c = —2 matter coupled to 2d gravity (to be 
described in detail in Sec. 4). Shown in red are also the shortest non-contractible 
loops, which we will discuss later. In addition Fig. 2 captures qualitatively the 
fractal structure of the ensemble of triangulations, to be analyzed below. 

The method described above for calculating the moduli of a genus one surface 
(and the moduli of a genus one triangulation) has a relatively simple generaliza- 
tion to higher genus surfaces and triangulations. This is outlined in appendix A. 
Below we describe the result of computer simulations which extract the moduli 
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of genus one triangulations. It is unproblematic to do the same for genus h trian- 
gulations, h > 1, although the moduh space now has the (real) dimension 6h — 6, 
but we have no continuum expressions with which we can compare the sampled 
distributions of moduli parameters, as already mentioned. 

3 Measuring r for c = 

In order to measure r for c = we perform Monte Carlo simulations starting out 
with a triangulation of the torus, constructed from N triangles. For a description 
of how to perform such simulations we refer to [14] . The Monte Carlo simulations 
consists of "moves" which change the triangulation locally, preserving N and the 
topology. After sufficiently many local moves we will have created a statistically 
independent triangulation. The number of moves needed on average for this 
increases with N. From the independent configurations constructed this way we 
now calculate r by 

(1) identifying two non-contractible link-loops 71 and 72 which generate the 
fundamental group of the torus, 

(2) constructing the harmonic differentials by solving the divergence and 
curl equations, = and {S(f))ijk — 0, 

(3) forming linear combinations such that (22) is satisfied, 

(4) constructing r using (13), 

(5) finding the modular transformation (11) which maps r to the fundamental 
domain A4 defined in eq. (12). 

In this way we construct a distribution of r's which we can compare with the 
theoretical distribution. The most computer intensive part of the calculation is 
step (2), the identification of the two-dimensional subspace of (discrete) harmonic 
1-forms. Fig. 3 shows a plot of the r-distribution. 

The fundamental domain A4 is non-compact. If we place the flat torus as 
shown in Fig. 4 we clearly have a torus very elongated in the vertical direction 
for large imaginary values of r. 

Of course this does not necessarily reflect a property of the original DT trian- 
gulation since we are using a (discrete) conformal map to embed the triangulation 
in the Euclidean plane. However, on average one would expect a large imaginary 
value of r to reflect a DT-torus with short non-contractible loops in one "di- 
rection" compared to the possible lengths of non-contractible loops in the other 
"direction". For the triangulations there are two restrictions entering: by defi- 
nition there is a shortest allowed loop length since we count the loop length in 
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Figure 3: (Part of) the fundamental domain. The left figure shows the density 
obtained from (a small number of) actual measurements for 1000 triangles, the 
right figure the theoretically calculated density. 

an integer number of links. Also, for a finite N, there is a limit to how long one 
can make geodesic loops. Thus we expect for a finite N that our DT-distribution 
of r's will fail for large imaginary values of r, but it should improve if (I) we 
increase and (II) we use DT-ensembles which allow smaller non-contractible 
loops. The second point can easily be tested in the following way: let us start 
with a DT-ensemble of "regular" triangulations DT{3). By that we mean trian- 
gulations where each vertex has at least order three and where two vertices are 
connected by at most one link. The dual graphs are 0^ graphs without tadpole 
and self-energy subgraphs. On this ensemble a shortest non-contractible loop is 
of length 3. We could instead allow two vertices to be connected by two different 
links, belonging to different triangles. We call this ensemble DT{2). The dual 
graphs are 0^ graphs without tadpoles, but with self-energy subgraphs (except 
the simplest ones which would result in vertices of order two). The shortest non- 
contractible loop length for this ensemble is clearly 2. Finally we consider the 
ensemble where we allow a link to loop to its own vertex. In terms of dual graphs 
we now allow tadpole graphs (but exclude still the simplest self-energy graphs 
to maintain that all vertices in the triangulation have an order larger than or 
equal three^). We denote this ensemble DT{1). The shortest loop length is one 

^We should stress that there is nothing important related to requirement that the vertex 
order is larger than or equal 3. We just maintained it because our original computer program 
had it built in. 
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Figure 4: The parallelogram rescaled and placed such that cui is the vector (1,0) 
and UJ2 = (ti, T2). 

for this DT-ensemble. In Fig. 5 we have shown the probabihty distribution of T2 
for a fixed N and compared it to the theoretical, continuum distribution. The 
DT{1) distribution is significantly closer to the continuum distribution for large 
r2 as expected. We also checked explicitly that the improvement was due to 
non-contractible loops of length 1 by calculating the distribution of T2 for the 
sub-ensemble of DT{1) which had no non-contractible loops of length 1. The 
results then agreed with the results for the DT(2) ensemble. In the rest of this 
article we will thus only use the ensemble DT{1), which we for simplicity just 
call the DT-cnscmble. 

Let us now turn to point (I) above, the N dependence. Fig. 6 shows that 
the T2 distribution agrees very well with the theoretical distribution for the large 
value of N used, and Fig. 7 shows the deviation from the theoretical distribution 
for different A^'s. 

Thus there seems to be a nice convergence to the correct distribution for 
T2. However, let us try to understand in more detail the relation between the 
geometry of the triangulations and the moduli parameter r. Due to (discrete) 
conformal invariance the relation cannot be very direct. We observed above that 
the minimal length of a non-contractible loop was important for the r distribution 
(for fixed A^). Let us classify a triangulation according to the length L of its 
shortest non-contractible loop. In the piecewise linear geometry of DT such a 
loop will be a geodesic curve in the DT sense. Thus we expect that its length 
L scales anomalously with respect to the area A'^ [9, 10]. We can check this 
explicitly by determining shortest closed loops for random triangulations in the 
DT-ensemble using the method described in appendix B (see also [13]). The 
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Figure 5: The cumulative probability of the various DT-ensembles plotted against 
T2. On the plot it is difficult to distinguish the curves for DT{2) and the DT{1) 
sub-ensemble with no non-contractible loops of length 1. The DT{1) curve is 
clearly closest to the theoretical curve. 



result for the expectation value (i^)Ar is shown in Fig. 8 and indeed we find that 

{L)n--N^'^\ 4 = 4, (23) 

where dh is the Hausdorff dimension of the DT-ensemble. One can view the 
triangulation to the left in Fig. 2 as an illustration of relation (23), the length of 
the shortest non-contractible loop shown in red is of the order iV^/'^'', not of the 
order -y/iV. Further implications of eq. (23) have been discussed in a recent paper 
[13]. 

From the continuum formulas (4) and (5) we expect the r distribution to fall 
off like e"'^'^^/^ for T2 much larger than 1. This is indeed in agreement with Fig. 
6. However, if we split the DT-ensemble in subsets according to L, we observe a 
uniwersa/ dependence in terms of the "dimensionless" variable LjN^I'^^: 

P{N,L,T2)^e-^^-f', /3 = 1 + 13.7 j . (24) 

The constants 13.7 and 1.85 are the result of a best ffi to the data (see Fig. 
9) and are not that precisely determined. The important point (apart from the 
universality) is that the main contribution for large T2 comes from small L/N^^^. 
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Distribution of x^ 




Figure 6: The T2 distribution for = 64000, compared to the theoretical distri- 
bution. 

Thus statistically there is a clear relation between small L's in the triangulation 
and large r2's in the embedded Euchdean plane. And "small" can be quantified 
as follows: since (i/j = 4 a typical linear extension of the triangulation, measured 
in geodesic distances, will be A^^/^, and "small" means small compared to this 
linear extension. Finally, let us emphasize that the existence of a small L in 
a triangulation does not imply that T2 is large. For instance both independent 
holonomy directions could have small non-contractible loops. As an extreme there 
could even be non-contractible loops of length one in both directions. In this case 
the embedding maps (like the ones shown in Fig. 1 and Fig. 2) to the Euclidean 
plane could be extreme, having only a few large triangles and the rest of the 
triangulation concentrated on a very little area. In such a case T2 would not be 
large, rather we would just have a nice illustration of the power of local conformal 
transformations. When A^ is not too large we indeed observe such situations since 
the smallest values of L are not that rare. 

4 Measuring r for c = — 2 

We want to test eq. (4) for c = —2. As mentioned c = — 2 is chosen because 
it is easier to address numerically than other matter systems coupled to two- 
dimensional Euclidean quantum gravity. The c = —2 matter coupled to quantum 
gravity has several realizations. A minimal (p, g)-conformal field theory, p < q, 
p, q co-prime integers larger than 1, has central charge c = 1 — 6{p — q)"^ /pq and. 
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coupled to quantum gravity, the string susceptibility for surfaces with spherical 
topology is 7o(c) = 1 — q/p. Formally a {p,q) = (1,2) model will have c = 
—2 and 70 = —1. It does not really belong to the minimal models and there 
is no Kac table associated with {p,q) = (1,2). Nevertheless one can find a 
fermionic conformal field theory which formally can be identified with a (1,2) 
model [15]. It can (partly) be viewed as a topological field theory. This is the 
reason the c = — 2 matter coupled to quantum gravity is often called topological 
2d gravity. However, there exists another formal representation of the c = —2 
matter system coupled to 2d Euclidean quantum gravity. It also leads to 70 = — 1 
after integrating over the matter fields, and in fact the value 70 = —1 was first 
calculated using this representation in the DT formalism [16]. 

Consider the partition function (1) for the bosonic string in d dimensions. 
Since the correspond to d Gaussian fields we can perform the integration over 
these fields and we will obtain 

J V^X, e-^«'^) ~ ( det(-A^)) (25) 

Here denotes the Laplace-Beltrami operator on the background geometry (jab 
and the prime signifies that the zero mode has been removed when calculating 
the determinant of A^. If we consider d as a formal parameter, d = —2 is special 
since we just get the determinant itself, as one would have for a suitable fermionic 
system, and since d is the central charge for positive integer, we formally have a 
system where c = — 2. The DT-formalism tells us that we should represent the 
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Figure 8: Log-log-plot of the expectation value (L) ^ versus A^. The fitted curve 
is given by {L)n = 0.471A^^/^ 

regularized partition function (1) as 

Z^'\^^) = J2^^''''^^^^^'Ti-A'^), (26) 

where the summation is over all triangulations of genus /i in a suitable DT- 
ensemble. /i denotes the bare cosmological constant of the 2d quantum gravity 
theory (the Liouville cosmological constant in the continuum notation), N{T) 
the number of triangles in T and Ct is symmetry factor of the triangulation T, 
i.e. the order of the automorphism group of the graph T. Finally At denotes 
the (discretely defined) Laplacian on the DT-surface, which we take to be the 
usual graph Laplacian of the (j)^ graph dual to the triangulation. We can trade 
the cosmological constant /i for the area of the (triangulated) surface. Since each 
triangle has the same area in the DT formalism, this implies keeping the number 
of triangles fixed and our partition function becomes 

zW(iV)= ^detT(-A^), (27) 

where DT{N) denotes the DT-ensemble with triangles. We expect Z^^\N) to 
approach the continuum Z^^\A) for iV — )■ oo via the identification A = Na^y/S/A, 
where a denotes the length of the links in the triangulation. 

It is well-known that the determinant of a graph Laplacian (with the zero- 
mode removed) is equal to the number of spanning trees on the graph. ^ Therefore 

■^For a given connected graph, a spanning tree is a connected subgraph which contains all 
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detT(— A^) = JV{T) where //{T) denotes the number spanning trees of the 0^ 
graph dual to the triangulation T. It was this representation of the combinatorial 
Laplacian which allowed to authors of [16] to solve the model for h = and prove 
that 7o = — 1, a result which was controversial at the time it was published. 
The partition function becomes 



where Tn{ST) denotes the set of genus h triangulations which are "decorated" by 
spanning trees. Thus two such triangulations are counted as different even if they 
as triangulations are identical, but if they have different "embedded" spanning 
trees. 

In [8] it was realized that this formula, in the case of genus zero {h = 0), could 
be used in computer simulations to generate directly, by a recursive algorithm, a 
set of graphs T which have the correct weight, including the matter fields. Thus, 
in order to test the properties of c = — 2 matter coupled to 2d quantum grav- 
ity, one can directly generate a set of independent graphs of arbitrary size and 
circumvent the problem of generating statistically independent configurations by 
local Monte Carlo updating, a problem which becomes increasingly time consum- 
ing when N is large. This way of dealing numerically with the c = —2 system 
has been used extensively to study the fractal properties of quantum gravity [17]. 

the vertices, but which contains no loops, i.e. closed paths of links. 





(28) 
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However, both the original solution of the model and the successive applications 
used in a crucial way that h = 0. For our application we need a generalization 
of the h = method to higher genera, or at least to genus one. Such a method 
is described for the simplest case h — 1 in the next section using a very recent 
result from [19]. In appendix A we sketch the procedure for larger genus. 

4.1 Random generation of decorated torus triangulations 

We will now describe an algorithm which generates a random triangulation T 
of genus h with triangles decorated with a spanning tree. The algorithm is 

designed such that any decorated triangulation is sampled with equal probability. 
The class of triangulations we are considering is the most general one, i.e. wc allow 
for triangles to be glued to themselves and two different triangles are allowed to 
be glued along more than one edge as long as the triangulation remains connected 
and the resulting piecewise linear manifold remains homeomorphic to a surface 
of genus h. 

Given a triangulation together with a spanning tree on its dual graph, we 
consider the set of edges of the triangulation which are not intersected by the 
spanning tree. This set of edges forms a graph consisting of + 1 links and 
containing 2h loops. If we cut open the triangulation along these links, we obtain 
a triangulation of the disc with a boundary consisting of A" + 2 links. This 
triangulated disc is completely characterized by the structure of a trivalent tree 
with A^ internal vertices (see the top part of Fig. 10). To get back to the original 
triangulation, the boundary edges of the disc have to be glued pairwise. This 
suggests that any decorated triangulation can be obtained by combining a tree 
and a pairwise gluing. 

Let us make this a bit more precise. We can view the A^ vertices of the span- 
ning tree as located at the center of the triangles. The order of the vertices (the 
number of links to which they belong) can be one, two or three. By construction 
they are located at triangles where two, one or no links have been cut open. Let 
us add two, one or no "external" links to the vertices. One can visualize this as 
the external lines still being located in the triangle to which the spanning tree 
vertex belongs, and "pointing" to the triangle links which have been cut open. In 
this way the spanning tree has been extended to a tree where the A" vertices have 
become internal vertices of order three and where N + 2 external links have been 
added. Wc mark one of the external links in order to have a rooted tree. In this 
way the tree becomes precisely of the form of a binary tree. The external links 
are now, by construction, in 1-to-l correspondence with the boundary links of the 
disc. The pairwise gluing of the edges of the disc to a genus h surface corresponds 
to a pairwise identification of the external lines in the binary tree such that the 
resulting trivalent graph can be placed on a genus h surface without any lines 
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Figure 10: The ingredients that specify a decorated genus triangulation with 
= 10 triangles. The top figure shows how a triangulated disc with a marked 
boundary edge corresponds to a rooted binary tree. The bottom figure shows 
how the pair-wise gluing a = {{1, 12}, {2, 11}, {3, 4}, {5, 6}, {7, 10}, {8, 9}} of 
the edges corresponds to a unicellular map which for genus is just a planar 
tree with a marked half-edge (i.e. the one labeled 1). To construct a genus 1 
unicellular map we select three vertices in the planar tree (the encircled ones). 
The distinguished half-edges are ai = 2, 02 = 3 and = 12. Therefore we should 
relabel {1,2,3,4,5,6,7,8,9,10,11,12} ^ {1,2,12,3,4,5,6,7,8,9,10,11} in a, 
yielding the new pair-wise gluing {{1, 11}, {2, 10}, {3, 12}, {4, 5}, {6, 9}, {7, 8}}. 

crossing (this is the meaning of the trivalent graph being of genus h). Such a 
pairwise gluing of a polygon is known as a unicellular map [19] (or one-face map) 
of genus h with + 2 "half-edges" . 

The above arguments show that we can generate a random decorated trian- 
gulation by separately generating a random trivalent planar tree and a random 
unicellular map. We then use the unicellular map to connect the external lines 
of the trivalent tree. 

As already remarked the random trivalent planar trees with a marked ex- 
ternal line are in 1-to-l correspondence with binary trees. There exists efficient 
algorithms to generate such trees, see e.g. [18] section 7.2.1.6. 

The major problem is to implement the random unicellar map. The c = — 2 
model was originally solved for /i = because it was understood that a genus zero 
unicellular map with A^ + 2 half-edges is given simply by a planar tree with one 
half-edge marked [16]. Fig. 10 illustrates this: the identification of the half-edges 
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has to form "rainbow" diagrams in order that one creates a genus zero surface. 
The corresponding planar trees are again related to binary trees and can thus be 
easily generated randomly. 

To generate a torus we need a random genus one unicellular map with N + 2 
half-edges. Luckily in [19] an explicit connection was found between unicellular 
maps of genus g and genus g + 1. In particular, for a genus zero unicellular map 
a procedure is given in which three distinct vertices are identified and the half- 
edges are relabeled in such a way that one obtains a genus one unicellular map. 
It is shown that any genus one unicellular map can be obtained through such a 
procedure in exactly two different ways (see [19], proposition 1 and corollary 1). 

Let us briefly summarize the procedure (sec Fig. 10 for an example). We 
label the half-edges of the N + 2-gon anti-clockwise by l,2,...,A^ + 2 and we 
provide them with the corresponding orientation, such that they have a starting 
vertex and a final vertex. A unicellular map of genus zero is fixed by giving a list 
a of A^/2 -I- 1 pairs of integers which tell us which edges to glue. After the gluing 
we have a tree with N/2 + 2 vertices (see Fig. 10) of which we randomly select 
three distinct ones. For each of them we select from the set of half-edges having 
that vertex as its final vertex the smallest index. We denote these indices by ai, 
a2, and as and reorder them such that l<ai<a2<a3<A'^-|-2. The resulting 
unicellular map of genus one is now given by gluing according to a in which we 
replace i — >■ f{i), where 



We refer the [19] for the actual proof of this statement. 
4.2 Numerical results for c = —2 

We use the above described algorithm to generate an ensemble of c = — 2 graphs 
and wc use this to perform the same measurements as for c = 0. Let us record the 
T2 distribution. It is shown in Fig. 11. We have a perfect fit even for a relatively 
small triangulation of 8000 triangles. In fact we have seen no deviation from the 
theoretical curve not compatible with the error-bars, so we have not addressed 
the approach to the theoretical curve as a function of N as we did for c — 0. 

We measured the expectation value (L) n of the length of the shortest non- 
contractible loop and the result is shown in Fig. 12. As for c = we expect 
anomalous scaling according to (23), only now dh is no longer equal to 4, but 
given by the formula [20] 




(29) 



4(c) = 2 
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Figure 11: The T2 distribution for = 8000, compared to the theoretical distri- 
bution (the red curve). 

for c < 1. The data shown in Fig. 12 is in perfect agreement with (30), which for 
c = — 2 becomes dh = (3 + a/T7)/2 f» 3.56. Eq. (30) was earher verified with good 
precision for genus zero surfaces by directly measuring the area N{r) enclosed 
within a circle of geodesic radius R and showing that {N{R))n ~ R'^'', with dh 
given by (30) [17]. As already noticed, (23) can be viewed as an independent ver- 
ification of the anomalous scaling of geodesic distance. Another rather stunning 
confirmation that the intrinsic structure of 2d quantum gravity is governed by the 
dimensionless quantity R/A^^'^'\ A being the area of the 2d universe, R a geodesic 
distance, can be obtained by looking closer at the actual probability distribution 
Pn{L) of the length of the shortest non-contractible loop as a function of N . In 
Fig. 13 we have shown the distribution for a range of N stretching from 126 to 
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Figure 12: Log- log-plot of the expectation value (L) ^ of the length of the shortest 
closed loop as a function of the volume N. The fitted curve corresponds to 
{L)j\i = 0.454A^-^/^'^^ (error-bars too small to display). 



more than 10^. In Fig. 14 we show that all of these Pn{L) are well described by 

P^iL) = N'/''^ Pix), x = ^. (31) 

This finite size scaling over such an amazing range of A^'s explains why we could 
hardly see any deviation from the continuum theoretical result in Fig. 1 1 : already 
= 8000 is in some sense very close to the N = oo limit according to (31). 




Figure 13: The distributions Pn{L) for = 126 up to 1024000 (the error bars 
are too small to display). 
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Figure 14: The rescaled distribution P{x) defined in (31) for = 126 up to 
1024000. 



For a conformal field tlieory witli central charge c we expect from the con- 
tinuum formulas (4) and (5) that the r distribution will fall off like e-'^i^-^)'^2/^ 
for T2 much larger than 1. We already checked this for c = and from Fig. 11 
it is clear that it is also true for c = —2. Again, if we split the DT-ensemble 
in subsets according to L, we observe a universal dependence in terms of the 
"dimensionless" variable L/N^^'^'^ (Fig. 15): 

P{N, L, T,) ~ e-f - ^ /3 = 1 + 11.4 (^^^ . (32) 

The constants 11.4 and 2.2 are different from the c = ones, and we have included 
a "shift" in the integer values L by a = 0.4, but the message is the same: at the 
boundary r2 = oo of moduli space the r distribution is completely determined by 
the smallest L's, where the independent statement of "small" is that L/N^/'^'^ 
is small. 

^The "shift" a can be viewed as a simple way to compensate for discretization effects for 
small L's. We expects a formula like (24) or (32) to reflect a continuum dependence Ij A^l'^^ 
where I is the continuum length t — La and A is the continuum area A (x Na^. While i can 
be arbitrarily small, this is of course not the case for L which is an integer. A priori there 
will be discretization effects if L is not much larger than 1. However, it is know that a shift 
L/N^/'^f^ {L — a)/N^^'^'^ can reduce the discretization effects [10, 21] and for c — —2 it does 
improve the fit. For c = the shift is not important and we left it out in (24). 
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Figure 15: The exponential fall-off P{N, L, T2) oc e '^^'^2/2 ^ _ _2^ formula 
(32). The fit corresponds to /3 = 1 + 11.4 (j0^) • 



5 Discussion 

Starting from a path integral approach to two-dimensional Euclidean quantum 
gravity coupled to matter fields, we need a regularization of the path integral if we 
want to have a definition of the theory which is not only formal. The formalism of 
dynamical triangulation is such a regularization. One can view the Z)T-formalism 
in two ways. In the first a triangulation represents a continuum piecewise linear 
geometry. Thus the geometry is considered fiat except at the vertices. At the 
vertices geometric quantities such as curvature, which involve derivatives of the 
geometry, can be singular. In a natural way the curvature has a delta-function 
distribution, being located at the vertices. One can take the point of view that the 
vertex represents a conical singularity with a certain deficit angle, but it might 
not be important to take such a literal continuum interpretation. At least it is 
a cumbersome road to take since we are interested in a limit where the number 
of vertices goes to infinity, and since we consider the triangles as building blocks, 
keeping them all as equilateral triangles, the deficit angles are not going to zero 
even if we consider triangulations with an increasing number of building blocks 
and even if we rescale all side lengths a of the triangles to zero. Alternatively we 
might simply consider a triangulation as a lattice realization of a 2d geometry, 
where not too much emphasis should be put on the piecewise linear structure. 

The analogy with the textbook derivation of the ordinary path integral in 
quantum mechanics might be useful: the path integral J Vx{t) ■ ■ ■ is obtained 
by discretizing the time interval, which leads to a multidimensional integral 





We can now 
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piecewise linear path x{t) connecting the initial point and the final point Xf. 
With such a choice, the piecewise linear paths with N vertices becomes a subset 
of the full set of continuous paths connecting Xi and Xf, which enter into the 
continuum path integral. It can be shown that this subset of piecewise linear 
paths is a dense set in the limit N ^ oo when the right norm is used to de- 
fine distances between continuum paths, namely the norm compatible with the 
Wiener measure (see [14] for a detailed discussion). However, we are not forced 
to take such a point of view for the particle path integral. We could simply 
view YliLi dx{ti) ■ ■ ■ as a lattice version of the formal continuum expression and 
whatever action we use in the continuum, we choose a suitable discretized lattice 
version of derivatives etc. Then, when the lattice spacing goes to zero we expect 
to obtain the continuum theory. Universality in the Wilsonian sense is then the 
key ingredient for obtaining a universal continuum limit. For the free relativistic 
particle this was analyzed in detail in [14] and the universality can be shown to 
be a consequence of the central limit theorem. 

Presently we have no rigorous mathematical definition of the formal measure 
J T^lg] over geometries but it is natural to believe that it involves the integration 
over all continuous geometries^. If we take the first point of view advocated for 
the path integral of the particle, we can choose to consider the DT-ensemble as a 
subset of the continuous geometries which, when the number of triangles N goes to 
infinity while the length a of the links goes to zero in such a way that the V — No? 
is kept fixed, hopefully becomes dense in the set of continuous geometries with 
volume V . However, contrary to the situation for the particle path integral, 
we cannot prove this since we do not presently know the norm which defines 
distances in the space of continuous two-dimensional geometries. As mentioned 
in the Introduction there is indeed some evidence that the DT-set of piecewise 
linear geometries can be viewed as dense in the set of continuous geometries, since 
a number of global quantities calculated using the formal continuum path integral 
agree with the corresponding quantities calculated using the DT-ensemble. 

However, the actual use of the £)T-ensemble, both in analytical calculations 
and numerical simulations, is more in the spirit of the second point of view ad- 
vocated for the path integral of the particle. According to this point of view we 
consider the DT-formalism as some conveniently chosen UV lattice representation 
of the two-dimensional geometries, where one should not put too much emphasis 
on the detailed piecewise linear geometric interpretation. Indeed, in the actual 
numerical simulations one measures geodesic distance using lattice links or dual 
links, not the actual geodesic distance in the piecewise linear geometry. Also, 

^Clearly, using the standard, formal, continuum procedure summarized in (1), one ends up 
with an integration over the field (p. The path integral over clearly involves fields (j) which 
are continuous but nowhere differentiable. The same will then be true for the corresponding 
geometries defined by the metric ga0 = e'^Qafi, where is a fixed background metric. 
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one is using discrete differences as substitutes for derivatives, rather than trying 
to formulate the field theory on a geometry with conical singularities. Similarly, 
one often considers discrete spin systems, like the Ising model, representing some 
conformal field theories. Such systems are of course only defined on a lattice. All 
numerical evidence as well as analytical calculations support the Wilsonian idea 
that lattice details are relatively unimportant when taking the continuum limit. 

It is with this second point of view in mind that we have analyzed the mod- 
uli structure provided to us by the DT-ensemble. We have not tried to view 
the triangulation as piecewise linear geometry with conical singularities to which 
one can associate a precise conformal structure as for instance described in [22] . 
Rather, we have just imitated a standard analytic construction of the moduli 
parameter r by replacing the harmonic differentials with suitable discrete differ- 
entials. This is in the spirit of the lattice approach. The conclusion: it works 
almost perfectly. In particular the c — —2 results, where we were not depend- 
ing on Markov chain Monte Carlo simulations, are in perfect agreement with the 
continuum results. Our numerical experiments also provide us with additional 
evidence that the anomalously scaling geodesic distance plays a key role in any 
deeper understanding of observables in quantum gravity theories. 
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A Generalization to higher genus 

The way we have set up the moduli measurements in section 2 allows for a rather 
straightforward generahzation to surfaces of genus h larger than 1.^ For such sur- 
faces the complex structure is characterized by the Teichmiiller space of complex 
dimension 3h — 3. Identifying explicitly 3h — 3 complex moduli parameters is a 
hard task and one which we will not pursue for general genus. Instead we will 
use a direct generalization of the genus 1 moduli parameter r in the complex up- 
per half- plane H to matrices in the /i-dimensional Siegel upper half- plane H^, 
known as period matrices [24]. The Siegel complex upper half- plane consists 

method similar to the one described here has been previously employed in a completely 
different context in [23]. 
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of all symmetric complex h by h matrices Q that have a positive-definite imag- 
inary part. A period matrix corresponding to a surface is known to completely 
determine its complex structure. For genus h < 2 we can actually identify 
with Teichmiiller space, but for /i > 3 Teichmiiller space appears as a non-trivial 
submanifold sitting in Mh (as is apparent from comparing their dimensions). 

In order to define the period matrix for a Riemannian manifold we have to 
first choose a set of 2h closed curves ai, . . . ,ah,bi, . . . ,bh which generate the 
fundamental group and which satisfy 

i{ai, Qj) — i{bi, bj) — and i{ai, bj) — 5ij (33) 

where denotes the oriented intersection number. The space of harmonic 

1-forms is 2/i-dimensional and it is possible to choose a basis {a^, Pi} dual to the 
curves Oj and 6j in the sense that 

[ Pj^ [ = and [ = S^j. (34) 

J ai J hi J ai Jbi 

A complex basis of holomorphic 1-forms is given by 

uji — ai + i (35) 

where * is the Hodge dual. The period matrix Q is then given by fl = A~^B in 
terms of the matrices 

^ and Bij = / a;,. (36) 

Jaj Jbj 

We can express these integrals in terms of the inner-products {ai\aj), {ai\f3j) and 
{f3i\f3j) by using the Riemann bilinear relations [24], which state that any two 
closed 1-forms p and a satisfy 

[ pAa = J2 [ P [ ^- [ ^ [ P- (37) 

J i=l 



Therefore 



Aij = Sij + i *Q;j = 5ij + i *ai A Pj = 5ij — i{ai\/3j) (38) 

Jaj J 



and 



Bij = i 



/ *ai = i / aj A *ai — i{ai\aj). (39) 

Jbi J 
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Another consequence of (37) is that 



J J j^—iJa^. Jbk J ak Jbk 

h 

= J2{ai\ak){Wj)-{ai\Pk){akm. (40) 



fc=i 



These formulae we can directly apply to DT piecewise linear geometries by 
replacing ctj and /5j by their discrete counterparts and the inner-product by the 
discrete one from (17). For genus h — 1 we used (40) to estabhsh the normaliza- 
tion of the discrete inner product leading to the expression for Q = [r] in (13). 
For genus h > 2 there is an overall ambiguity in the definition of Q because 
(40) is not exactly satisfied (up to overall factor) due to discretization artefacts. 
However, we expect (and we have confirmed this numerically for genus 2) that 
for large random surfaces (40) is close to a multiple of the identity matrix with 
high probability. In that case we can unambiguously normalize the inner-product 
and determine Q. 

The modular group SL{2,Z)/Z2 which acts on the upper half-plane as in 
(11) generalizes to the action of the symplectic group Sp{2h, Z)/Z2 on the Siegel 
upper half- plane Mh- Fundamental domains can again be worked out but become 
increasingly cumbersome for larger genus (see e.g. [25]). 

Finally let us mention that our algorithm for generating random DT surfaces 
coupled to c = —2 conformal matter can be straightforwardly extended to genus 
h >2. The only missing ingredient is a construction of a random unicellular map 
of genus h. Here again we can use results from [19] in which for any genus h 
an explicit bijcction is found between the set of unicellular maps of genus h and 
a union of sets of unicellular maps of lower genus with a particular number of 
distinguished vertices. 



B Finding shortest non-contract ible loops 

First let us present a method of constructing curves in the triangulation that 
generate the fundamental group. Inspired by the methods described in Sec. 4.1 
we generate an arbitrary spanning tree on the (p^ graph dual to the triangulation. 
Then we consider the set of edges of the triangulation that are not intersected 
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by this spanning tree/ As mentioned in Sec. 4.1 the graph G formed by these 
edges contains 2h cycles, where h is the genus of the triangulation. To extract 
these cycles we choose an arbitrary vertex v m G and generate a spanning tree 
for G based at v. This tree will contain all but 2h of the edges of G. Adding 
any of the remaining 2h edges to the tree will lead to a cycle and therefore to a 
unique closed path based at v. By rearranging (and if necessary composing) the 
paths thus obtained we arrive at a canonical set of generators {7j}j=i,...,2/i, i-e. a 
set that satisfies (33), for the fundamental group of the triangulation. 

Once we have such a set of generators we can establish for any closed curve 
whether it is contractible or not by computing its oriented intersection number 
with the generators 7j. The curve is contractible if and only if all these intersec- 
tion numbers vanish. To make this test more efficient notice that we can easily 
construct by hand for each generator 7^ a closed discrete 1-form 0j, i.e. satisfying 
d4>i = 0, such that the intersection number of a curve with is equal to the 
discrete integral of (pi along that curve. Hence, a closed curve is non-contractible 
if the integral of at least one of the 1-forms 0i is non- vanishing. 

Given a vertex v we can find a shortest non-contractible loop based at v by 
performing a so-called breadth-first search in the edge-graph of the triangulation 
starting at v. Once we encounter a vertex that we have already visited before, 
we have implicitly established a loop in the edge-graph. The first such loop we 
meet that is non-contractible will automatically have minimal length. 

Now in principle we can repeat this procedure for each vertex v in the tri- 
angulation to find the overall shortest non-contractible loop (or rather a non- 
contractible loop of minimal length as there usually more than one) . However in 
general the set V of vertices for which we have to perform this procedure can be 
greatly reduced. Indeed, we known that any non-contractible loop will intersect 
at least one of the generators 7^, so it suffices to take V to consist of all vertices 
contained in the 7^. In order to obtain such a set V with as few vertices as pos- 
sible it is worthwhile to first spend some time to shorten the 7^. This will result 
in a set V with a number of vertices of the order N^^^>^ with N the number of 
triangles. Since a single breadth-first search involves a number of steps of the 
order A^, the full algorithm will have an expected run-time of the order N^'^^^^'^, 
which amounts to A^^-^^ for c = and A^^-^^^ for c = —2. 
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